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Q H i The kinematic parameters of the Gould Belt are redetermined using modern data 

q . on the motion of nearby young (logt < 7.91) open clusters, OB-associations and 

moving stellar groups. The modelling carried out shows that the residual veloci- 
ties achieve the maximal values equal to —4 km s _1 for rotation (in the direction 
of the Galactic rotation) and +4 km s _1 for expansion, when the distance from 
the kinematic centre is about 300 pc. We assumed the following parameters of 
the Gould Belt centre: l = 128° and R = 150 pc. It is shown that the whole 
structure of the Gould Belt is moving relative to the local standard of rest with a 
velocity of 10.7 ± 0.7 km s _1 in the direction towards I = 274 ± 4°, b = -1 ± 3°. 
Using the rotation velocity, we obtained the virial estimation of the Gould Belt 
mass as 1.5 x 1O 6 M . 
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The Gould Belt is a system of nearby OB-associations [1], HI clouds [2], H 2 
molecular cloud complexes [3]. At present, there exist several scenarios for 
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the formation of the Gould Belt. According to the first scenario, the Gould 
Belt is formed as a result of the collision of high-speed neutral hydrogen 
clouds with the Galactic disk [4]. According to the second scenario, the Gould 
Belt is formed owing to explosion of supernovae [2, 5]. According to the third 
scenario [6], the formation of the Gould Belt is a stage of kinematic evolution 
of the Local Stellar System. Linblad [7] suggested a model of proper rotation 
and expansion of the Gould Belt which takes into account the inclination of 
the disk to the Galactic plane (20°). In the work of Bobylev [8], the well- 
known Linblad approach is developed for the case when star heliocentric 
distances are determined exactly. 

The goal of this work is redetermination of kinematic parameters of the 
Gould Belt structure using modern data on open clusters. The advantage 
of such an approach is that the estimations of ages and distances of open 
clusters are more reliable than those of individual stars. 

The catalogue of 652 open clusters hereafter we denote as COCD [9,10]. 
Our working list is composed primarily from the COCD young open clusters 
with age logt < 7.9, denoted in the paper of Piskunov et al. [10] as OCC1. 
In the COCD catalogue, random errors of mean proper motions of the clus- 
ters are small because of the use of the combined catalogue ASCC-2.5 [9]. 
The radial velocities of the COCD open clusters are based on the catalogue 
published by Barbier-Brossat and Figon [11]. We use more comprehensive 
star radial velocities from the Orion Spiral Arm Catalogue (OSACA) [12], 
which combines a number of large catalogues of radial velocities including 
those of Duflot et al. [13], Barbier-Brossat and Figon [11], Nordstrom et 
al. [14] and Famay et al. [15]. In comparison with the work in [3], we made 
essential supplements for stars of the Sco-Cen complex using the list com- 
piled by de Zeeuw et al. [1] and the OSACA. We considered the open cluster 
Chamaeleontis. We used modern data for 25 stars of cluster TWA as well. 
Our approach to kinematic analysis is based on the following equations [8]: 

V r = -« cos b cos I — v © cos b sin I — w & sin b + cos 2 bk r + (1) 
+ (R — R Q )(r cos b — R cos(/ — l )) cosbk' — R Q (R — R ) sin(7 — l Q ) cosbuj' , 
4.7Arfii cos b = u Q sin / — v@ cos/ — (2) 
— (R — R )(RoCOs(l — l ) — rcosb)cu' + rcosbcu + R (R — R ) sin(Z — l ))k' , 
4.74r/^ = u Q cos / sin b + v Q sin / sin b — w Q cos b — cos b sin bk r — (3) 
— (R — R Q )(rcosb — R Q cos(/ — l )) sinbk' + R Q (R — R ) sin(/ — l ) sin bui' Q . 

Here r = 1/ir is the heliocentric distance of a star, R Q is the distance from 
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Figure 1. Residual velocities Vr (top) and Vg (bottom) versus distance R from 
the kinematic centre l Q = 128°, R Q = 150 pc; the vertical line indicates the value 
of R . 




the Sun to kinematic centre, l is the direction to the centre, u Q ,v Q ,w Q are 
the components of the solar peculiar velocity. The components cos b and 
fib are expressed in milliarcseconds per year (mas yr" 1 ), the radial velocity V r 
in kilometers per second (km s -1 ), the parallax n in milliarcseconds and the 
distances R, R Q and r in kiloparsecs (kpc). The quantity u a is the angular 
velocity of rotation, k Q is the radial expansion/contraction velocity of a star 
system at the distance R Q and parameters u' Q and k' are the corresponding 
derivatives. The distance R from a star to the centre is evaluated as R? = 
r 2 cos 2 b — 2R r cosbcos(l — l ) + R 2 ,. The system of conditional equations (1)- 
(3) contains seven unknowns which we determine using the minimum least- 
squares technique. It is assumed that the left-hand sides of the equations are 
free from the Galactic differential rotation. 

The second method of kinematic analysis is based on linear residual ve- 
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Figure 2. Smoothed distribution of [TV-residual velocities. 

locities which we approximate as: 

V l = V , i + a l (R-R ) + b l (R-R ) 2 , (4) 

where % = (R, 6), and the residual velocity Vr of expansion and the residual 
velocity Vg of rotation must be obtained by decomposition of the spiral ve- 
locity components U and V into radial and tangential parts using the known 
centre parameters l and R . The results of the analysis are presented in 
figures 1 and 2. We can conclude that it is shown in our work that the Gould 
Belt structure takes part in several motions. Firstly, the whole complex, 
free from common Galactic rotation, is moving relative to the local standard 
of rest [16] with a velocity of 10.7 ±0.7 km s _1 in the direction towards 
I = 274 ± 4°, b = — 1 ± 3°. Secondly, residual rotation and expansion of the 
whole system exist. We used the following parameters of the kinematic cen- 
tre: l = 128° and R = 150 pc [8]. The modelling carried out shows that the 
residual velocities are reliable in a small region near R « R Q and are equal 
to —2.8 ± 0.7 km s -1 for rotation (in the direction of the Galactic rotation), 
and to +3.3 ±0.7 km s _1 for expansion. The maximal values are achieved 
when the distance from the kinematic centre is about 300 pc and are equal 
to —4.3 ± 1.9 km/c for rotation and to ±4.1 ± 1.4 km s" 1 for expansion. 
Using the rotation velocity, we obtained the virial estimation of the Gould 
Belt mass as 1.5 x 1O 6 M . 

Acknowledgement 



4 



This work is supported by the Russian Fund for Basic Research (grant No 05- 
02-17047). 

References 

[1] P. T. de Zeeuw, R. Hoogerwerf, J. H. J. de Bruijne et al, Astron. J. 117 
354 (1999). 

[2] W. G. L. Poppel and P. Marronetti, Astron. Astrophys., 358 299 
(2000). 

[3] C. A. Perrot and I. A. Grenier, Astron. Astrophys. 404 519 (2003). 

[4] F. Comeron and J. Torra, Astron. Astrophys. 281 35 (1994). 

[5] C. A. Olano, Astron. Astrophys. 112 195 (1982). 

[6] C. A. Olano, Astron. Astrophys. 121 295 (2001). 

[7] P. O. Lindblad, Astron. Astrophys. 363 154 (2000). 

[8] V. V. Bobylev, Pis'ma Astron. Zh. 30 861 (2004). 

[9] N. V. Kharchenko, Kinemat. Phys. Neb. Tel, 17 409 (2001). 

[10] A. E. Piskunov, N. V. Kharchenko, S. Roser et al, Astron. Astrophys. 
445 545 (2006). 

[11] M. Barbier-Brossat and P. Figon, Astron. Astrophys. Suppl. Ser. 142 
217 (2000). 

[12] V. V. Bobylev, G. A. Gontcharov, and A. T. Bajkova, Astron. Zh. 83 
821 (2006). 

[13] M. Duflot, P. Figon and N. Meyssonier, Astron. Astrophys. Suppl. Ser. 
114 269 (1995). 

[14] B. Nordstrom, M. Mayor, J. Andersen et al, Astron. Astrophys. 418 
989 (2004). 



5 



[15] B. Famaey, A. Jorissen, X. Luri et al, Astron. Astrophys. 430 165 
(2005). 

[16] W. Dehnen and J. J. Binney, MNRAS 298 387 (1998). 



6 



